Quasispherical gravitational collapse in 5D Einstein-Gauss-Bonnet gravity 
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We obtain a general five-dimensional quasispherical collapsing solutions of irrotational dust in 
Einstein gravity with the Gauss-Bonnet combination of quadratic curvature terms. These solutions 
are generalization, to Einstein-Gauss-Bonnet gravity, of the five-dimensional quasispherical Szkeres 
like collapsing solutions in general relativity. It is found that the collapse proceed in the same 
way as in the analogous spherical collapse, i.e., there exists a regular initial data such that the 
collapse proceed to form naked singularities violating cosmic censorship conjecture. The effect of 
Gauss-Bonnet quadratic curvature terms on the formation and locations of the apparent horizon is 
deduced. 
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I. INTRODUCTION 

One of the remarkable special classes of the exact so- 
lution of the nonvacuum Einstein equations G a b = KT a j 
available to date was discovered by Szekeres [l[ . They are 
obtained by solving Einstein equations with irrotational 
dust T a b — eu a Ub as the source, for the line element 



ds 2 = -dt 2 + X 2 dr 2 + Y 2 (dx 2 + dy 2 ), 



(1) 



relative to comoving coordinates. Here u a is a veloc- 
ity (i.e. unit timelike) vector field and e is the energy 
density of the system. One can say that Szekeres so- 
lutions are obtained when the spherical symmetry orbits 
in the Lemaitre-Tolman-Bondi (LTB) model Q are made 
nonconcentric to destroy the symmetry, but the energy- 
momentum tensor is still that of dust. Therefore, the 
Szekeres spacetime is often called as quasispherical, for 
dcfmiteness we shall name it as quasispherical Szekeres 
(QSZ) solutions. Being an exact model of the spacetime 
geometry, the QSZ solutions have primarily been adapted 
with regards to studies of nonspherical collapse of inho- 
mogeneous dust cloud P, [3rJ §| in four dimension s (4D), 
and in higher dimensions [3, Qi| , in cosmology [Il1 - fl6| 
and also in observational cosmology [l7|- They are very 
important because they admit no Killing vectors and 
hence they are lacking symmetry. 

In any attempt to perturbatively quantize gravity as 
a field theory, higher-derivative interactions must be in- 
cluded in the action. Such terms also arise in the effective 
low-energy action of string theories. Among the higher 
curvature gravities, the most extensively studied the- 
ory is the so-called Einstein-Gauss-Bonnet (EGB) gravity 
[la - t29| . The EGB gravity is a special case of Lovelocks' 
theory of gravitation, whose Lagrangian contains just the 
first three terms. Gauss-Bonnet gravity provides one of 
the most promising frameworks to study curvature cor- 
rections to the Einstein action in supersymmetric string 



theories, while avoiding ghosts and keeping second order 
field equations. 

In Einstein-Gauss-Bonnet gravity less number of exact 
solutions have been known so far. Static and spherically 
symmetric black hole solutions with or without a cosmo- 
logical constant, as well as a topological black hole in an 
anti-de-Sitter spacetime were obtained [l8l - [22| [. The ef- 
fects of Gauss-Bonnet terms on the Vaidya solutions have 
been investigated in [23l - r27j . and on the LTB solutions 
in [2^, [29| . Here, we consider the five-dimensional (5D) 
action with the Gauss-Bonnet terms for gravity and give 
an exact model of the quasispherical gravitational col- 
lapse. Using our new solution, we investigate the nature 
of singularities of such a spacetime in terms of its being 
hidden within a black hole, or whether it would be vis- 
ible to outside observers and compare it with analogous 
relativistic case. Thus, the aim of this paper is to extend 
the previous studies on the quasispherical gravitational 
collapse of inhomogeneous dust, including the second or- 
der perturbative effects of quantum gravity solutions: a 
Gauss-Bonnet generalization of the QSZ solutions in 5D 
spacetime, namely, 5D-QSZ-EGB. 

There are several issues that motivate our analysis: 
how does the Gauss-Bonnet term affect the final fate of 
collapse? What is the horizon structure in the presence 
of the second order perturbative effects of quantum grav- 
ity? Whether such solutions lead to naked singularities? 
Do they get covered due to departure from spherically 
symmetry? Does the nature of the singularity changes in 
a more fundamental theory preserving cosmic censorship 

El? 



II. SZEKERES SOLUTIONS IN 5D EINSTEIN 
GAUSS-BONNET GRAVITY 

In this section, we derive the relevant equations for 
Szekeres model for irrotational dust to the 5D Gauss- 
Bonnet extended Einstein equation - 5D-QSZ-EGB so- 
lutions. We begin with the following 5D action: 
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where R is 5D the Ricci scalar, and K5 = y/SnG^ is 5D 
gravitational constant which is set to unity. The Gauss- 
Bonnet Lagrangian is of the form 



Lqb — R 2 ~ 4:R a bR a + R a bcdR a 



bed 



(3) 



where a is the coupling constant of the Gauss-Bonnet 
terms. This type of action is derived in the low-energy 
limit of heterotic superstring theory [30]. In that case, 
a is regarded as the inverse string tension and positive 
definite, and we consider only the case with a > in this 
paper. 

The action ([2]) leads to following set of field equations: 



Gab = G a b + aH a 



T, 



a b ; 



where 



Gab — Rab — ydabR, 



is the Einstein tensor and 



Hab — 2[RRab ~ 2R aa Rb 
^-jSab^GB, 



2R ap R, 



R 



a/37 



(4) 
(5) 
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is the Lanczos tensor. 

To find QSZ spacetime, Szekers [l| imposed no a priori 
symmetry assumption, but instead postulated the metric 
in a special form which, for our 5D case, can be written 
as: 



ds 1 



-dt 2 + X 2 dr 2 +Y 2 {dx 2 + dy 2 + dz 2 ). (7) 



The geodesic fluid flow vector is d/dt and the coordinates 
(x, y, z) are comoving spatial coordinate constant along 
each world-lines so that u a — 8f. X and Y are func- 
tions of (t, x, y, z) to be determined from the Einstein 
equations. The energy-momentum tensor for dust is 



Tab = e(t,x, y, z)5 t a S t b 



(8) 



It is seen that there exist solutions according as Y' = 
or as Y' ^ (throughout this paper . = d/dt and 
' = d/dr). In 4D, the family of the solution correspond- 
ing to Y' = is a coincidental generalization of the Fried- 
man and Kantowski-Sachs models. Bolejko et ah [l6| 
pointed out, in 4D, that the case Y' = found no useful 
application in astrophysical cosmology. Hence, we shall 
confine our discussion to the case Y' 7^ 0. After the Ein- 
stein equations are solved, it is required that there must 
exists two functions R(t, r) and u(r, x, y, z) such that 



Y = 



R(t,r) 



P(r,x,y,z) 



X 



PY' 
W(r)' 



(9) 



(10) 



Here, W = W(r) is an arbitrary function of r with re- 
striction that W(r) > 0. The solution for P reads 

P = A(r)(x 2 +y 2 + z 2 ) 

+B x (r) 2 x + B y (r) 2 y + B z {r) 2 z + C(r), (11) 

with the free functions A(r), B x (r), B y (r), B z (r) and 
C(r) satisfying an algebraic equation 



AC - Bi 



Bt 



Bi = -, 

z 4' 



e = 0, ±1. 



(12) 



The factor e determines whether the 2-surface is spherical 
(e = +1 ), p seudospherical [hyperbolic] (e = -1) or planar 
(e = 0) [lj]. Upon introducing transformations 



x = sin if) sin <j) cot — , 



y = cos ip sin (f> cot • 



cos cot 



2' 



(13) 



the 3-metric yields the more familiar spherical form 

R 2 {t,r)(d6 2 + sin 2 9(d(l> 2 +sin 2 0dV 2 ))- (14) 

As can be seen, if t = const, and r= const, the above met- 
ric becomes the metric of the three-dimensional sphere. 
Hence, every t = const, and r= const slice of the Szekeres 
spacetime is a sphere of radius R. Thus the term quasi- 
spherical adapted. However, the spheres under consider- 
ation are not concentric. Thus, the QSZ solutions con- 
sidered here are a generalization of the 5D-LTB solutions 
in which sphere of constant mass are made nonconcen- 
tric. The A(r), B x (r), B y {r), B z (r) and C(r) determine 
how the center of a 3-sphere changes its position in a t — 
constant space when the radius of the sphere is changed. 
When the spheres are concentric, the metric (T7|) becomes 
the line element of the 5D-LTB model [3ll - l34j which 
arises when we set A — C — 1/2, B x = B y = B z = 0, 
with no loss of generality. 

The acceleration equation is given by the Q r r = com- 
ponent of the field equations 



R 
R 



1 * (W 2 -l) „ R 2 



(W 2 - 1) 
~R 2 



R 2 



= 0. 



(15) 



Thus, it is necessary that R satisfies the Friedmann-like 
equation 



R- 



1 - 4a- 



W 2 



1 



R 2 



(W 2 -l) + ^-2a^. (16) 



Here, F — F(r) is an arbitrary function of r and is re- 
ferred to as mass function. In a Newtonian limit F(r) is 
equal to the mass inside the shell of radial coordinate r, 
assuming that the mass function has no angular depen- 
dence. The function F must be positive, because F < 
implies the existence of negative mass. Equation (|16p is 
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the master equation of the system which governs the dy- 
namical properties of the system, which is the same as 
one obtains in a spherically symmetric collapse of dust in 
EGB HI. Hence, we can refer to the the 5D-QSZ-EGB 
discussed here as quasispherical collapsing space-times. 

It is now straight forward to calculate the energy den- 
sity e, which we obtain by substituting the solutions for 
X and Y . The mass density is 



3 PF' - AFP' _ 3 PF' - AFP' 
e{t,x,y,z) - - p3xy3 - 2B?(PR' -RP'y 



(17) 



ical symmetry these solution corresponds to more con- 
versant 5D Tolman-Bondi models [32], [33[ . 

In many instance, it may be preferential to use in our 
work not density e, but mass M. M is mass within sphere 
(r=constant) at given time t, which we compute as fol- 
lows 



M(t,r) = J dr I I / dx dy dz XY 6 e. (IS) 



This solution has in general no symmetry. For the spher- Upon using the density Eq. (fl"7| in (fT8| , we get 



M(t,r) 



dr — 
2 



EL 
w 



PF'-AFP' 

dx dy dz 

y WP 3 

dx dy dz 2F d 

P 2 + ~W dr 



dx dy dz 
P 2 



(19) 
(20) 



Now dx dy dz / P 2 is a metric on the unit 3-sphere, so 



dx dy dz , 
% = 2tt 2 

p2 



so that 



and 



dM 
dr 



M(t,r) 



= 2tt z XY a e 



3vr 2 / — dr 
Jo W 



(21) 



(22) 



(23) 



This is a generalization of a formula obtained for mass 
function in 4D QSZ [1]. In addition, e > implies 
dM(r)/dr > or F'{r)W(r) > 0. 

It is easy to see that as a — > the master solution 
(UHJ) of the system reduces to the corresponding 5D-QSZ 
solution in 0, [Toj 



R 2 = W 2 



F 



(24) 



Rotation and acceleration of the dust source are zero, the 
expansion is 



h XL ,'1. 

Y> R 



and the shear is 
S 2 = 



•XL 



12 RR +9 & 



(25) 



(26) 



The 5D-LTB limits of these scalars can be obtained by 
replacing Y by R in the above expressions. 



We assume, as in QSZ models Q, the following regu- 
larity conditions:- 1. The metric is C . 2. There are no 
shell crossing singularities (Y' = PR'-RP' > 0). 3. The 
metric is locally Euclidean, i.e., we must have W(0) = 1. 

The dynamical properties of the model discussed above 
is governed by Eq. (fTB|) . which is the same as one gets 
in the corresponding spherical symmetric collapsing solu- 
tions [2^ |. To get the 5D homogeneous dust limit, we put 
W 2 (r) = 1 + k r 2 , F(r) = Kr 3 , where k = 0,±1 and 
K is constant. Then one gets three type of solutions - 
parabolic, elliptic or hyperbolic according as W 2 (r) = 1, 
W 2 (r) < 1 or W 2 (r) > 1. 

In particular, the 5D Friedmann-Robertson- Walker 
metric follows when R(t, r) — S(r)f(t) and W(r) = 
1 — koS 2 (r), where ko is the curvature index of 5D 
Friedmann-Robertson- Walker models and B x = B y = 
B z — 0, C = 4 A = 1. Then the metric obtains usual 
form 



ds 2 



-dt 2 



1 - k S 2 



dS 2 + S 2 f 2 (t)dQ 2 3 . (27) 



A. Solutions for the zero-energy case, W(r) = 1 

It may be noted that in the general relativistic case 
(a — > 0), Eq. (fT6|) has three types of solutions and 
W(r) > 1, W(r) = 1 or W(r) < 1 determines the type 
of evolution. From Eq. (I16|) , we obtain 



R 2 = (W 2 - 1) 



1T\ A 



16a 2 



(w 2 



8aF(r) 
i? 4 



(28) 
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There are two families of solutions that correspond to 
the sign in front of the square root in Eq. (|2"51) . We 
call the family that has the minus (plus) sign the minus 
(plus) branch solution. In the general relativistic limit 
a — > 0, we recover the 5D-QSZ solution in Einstein grav- 
ity @, Maeda [H has analyzed LTB models near 
center (r ~ 0) in EGB and pointed out the occurrence 
of major changes in the final fate of collapse (see also, 
_28]). Here we present the 5D-QSZ-EGB exact solution 
in close form, which facilitates us to analyze the final 
fate of gravitational collapse. The condition W{r) = 1, 
is the marginally bound condition, meaning the collaps- 
ing shell is at rest at spatial infinity (R — oo) with zero 
energy in the infinite past. In the present discussion, we 
are concerned with gravitational collapse, which requires 
R(t,r) < 0, i.e., we assume that all portion of the dust 
cloud are momentarily collapsing. Eq. (|28| can be inte- 
grated to 



2y/2 



tan 



3R 2 - VR 4 + 8aF 
2V2R[VR T TS^F - R 2 } 1 / 2 



aR 2 



VR* + SaF - R 2 



(29) 



where t^(r)(r) is an arbitrary function of integration. 
Here we note that the solutions R(t, r) are same as the 
5D-LTB-EGB models [H[ and are not affected by the 
dependence of the Szekeres model on the (x, y, z) coor- 
dinates. As it is possible to make an arbitrary relabel- 
ing of spherical dust shells by r — > g(r), without loss of 
generality, we fix the labeling by requiring that, on the 
hypersurface t = 0, r coincide with the area radius 

R(0,r)=r. (30) 

This corresponds to the following choice of t q (r) : 



us assume at t = t c (r), we have R(t, r) = 0, which is the 
time when the matter shell r — constant hits the physical 
singularity, i.e., t c (r) is defined by R(t c (r),r) = 0. 

The central singularity curve can be obtained using 
Eq. O as 



(33) 



which represents the proper time for the complete col- 
lapse of a shell with coordinate r. Interestingly, posi- 
tive a delays the formation of singularity. In the limit 
of vanishing a we recover the crunch time for relativistic 
5D-QSZ-EGB. The eight arbitrary functions A(r), B x {r), 
B y {r), B z (r), C(r) F(r), W(r) > and t c (r) completely 
specify the dynamics of collapsing shells. However, only 
seven are independent because relation (|12[) is between 
them. With the choice of r, one can fix one more func- 
tion and thus the total degree of freedom is six. The 
corresponding 5D-LTB-EGB models have only three free 
functions: F(r), W(r) and t c (r). Hence the 5D-QSZ- 
EGB is functionally more generic. 

In order to study the collapse of a finite spherical body 
in EGB, we have to match the 5D-QSZ-EGB solution 
along the timelike surface at some r = r c > to the 5D- 
EGB Schwarzschild exterior discovered by Boulware and 
Desser [l8[, and Wheeler [2(|. The analysis is similar 
to the case of matching5D-LTB-EGB to the 5D-EGB 
Schwarzschild exterior [29J and, will not be discussed 
here. Bonnor [ll| was the first to proved that Szek- 
eres spacetime can be matched to Schwarzschild vacuum 
spacetime. 



III. APPARENT HORIZON AND TRAPPED 
SURFACE 



2V2 



tan 



3- Vl + SaF 
2V2[ a/I + 8aF - 1]V2 



Vl + SaF - 1 ' 



(31) 



where F = F/r A . Now, we discuss the nature and oc- 
currence of curvature singularities in the 5D-QSZ-EGB 
models. It follows from Eq. (fTT)) that the curvature sin- 
gularities occurs in 5D-QSZ-EGB solutions when 



R = or Y' = PR 1 - RP' = 0. 



(32) 



Szekeres [l[ analyzed the singularities of the first kind 
and second kinds according to the above two conditions. 
As is the case of the LTB, in 4D, Y' = corresponds to 
a shell crossing, however it is qualitati vely different from 
that which occurs in LTB models ( see [la] for details) . In 
this paper we shall restrict ourselves to the singularities 
of the first kind which we call the central singularity. Let 



The apparent horizon (AH) is the outermost 
marginally trapped surface for the outgoing pho- 
tons. The AH can be either null or spacelike, that is, it 
can "move" causally or acausally. The main advantage 
of working with the apparent horizon is that it is 
local in time and can be located at a given spacelike 
hypersurface. The event horizon instead is nonlocal. 
Moreover, the event horizon is a an outer covering 
surface of apparent horizon and they coincide in case 
of static or stationary spacetime. Trapped surface is 
defined as a compact spacelike 2-surface both of whose 
future pointing null geodesies families are converging. 
Physically, it captures the notion of trapping by imply- 
ing that if 2-surface S7 r ,t) (t, r = contant) is a trapped 
surface then its entire future development lies behind a 
horizon. To obtain criterion for existence of such S< r ,t) 1 
let K a be the tangent vector to the null geodesies. It 
follows that along null geodesies, we have 



K^K» = 0, 



K? V K V = 0, 



(34) 
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and K x = K y = K z = 0. For the metric the above 
defined null congruence satisfies the following condition: 

{K 1 ) 2 - X 2 {K r ) 2 = 0, (35) 

on S( r> t)- A choice of affine parameter 

K* = X, K r = e, (36) 

with e — ±1, clearly satisfies the condition (f3"5"j). The 
positive (negative) sign of invariant Kt^ is determines 
the divergence (convergence) of the null geodesic. Also 
± corresponds to the outward and inward geodesic re- 
spectively. Now, 





Y'\ 























K^ = K% + K^ r 

Upon taking the time derivative of the null condition 
K^Kf, = 0, we get 



K\ - eXK\ - X = 0. 



(38) 



In the second condition Ktl K u — 0, consider the // = 1 
component; 



VI 

K r r + eXK\ + e— + 2X = 



(39) 



Eliminating K r t from the previous two equations we can 
rewrite Eq. (|37|) as 



Y> y 



(40) 



Thus the apparent horizon must satisfy 

eY' +17 = 0. (41) 

In absence of shell crossings (Y' ^ 0), and using Eqs. © 
and (pTJ]) . we have 



R 2 {t AH {r),r) =-W(r). 



(42) 



Considering Eq. (|16j) . the apparent horizon condition 
(g2|) becomes 



fl(t AH (r),r) = y/F(r) - 2a. 



(43) 



One can also employ the usual condition for the existence 
of the apparent horizon 



9 ab Y a Y b = 0. 



(44) 



Considering Eqs. (f9]) . ([TUT) and (fT6| . the apparent horizon 
condition (|44|) again gives the same result (|43|) . In the 
relativistic limit, a — > 0, Rah — > \fF{ r ) ,32]. Further, 
Eq. (j43|) has a mathematical similarity for the analogous 
situation in null fluid collapse where the expression for 
apparent horizon is tah = \J m(v) — 2a [27| . 



IV. END STATE OF COLLAPSE 

In this section, we analyze end state of the collapse of 
5D-QSZ-EGB dust collapse in terms of the given regular 
initial density and velocity profile. We denote p(r) as 
the initial density of the dust cloud at t = for a fixed 
radial direction (x — y — z — constant). It is assumed 
that p, P and F to be expandable around r = on the 
initial hypersurface. Further the function P is restricted 
by Eqs. (|TT|) and (Tj"2"]) . Henceforth, in this section, we 
adopt here a method similar to 0, Q which we modify 
here to accommodate the higher dimension spacetime. 
We take 



n=Q 



P = J2 P nr n , 



(45) 



(46) 



71=0 



where Pn ^ and the regularity condition implies that 
Pi = [|, and 



F = J2 F nr n 



+ 4 



(47) 



n=0 



Therefore, 

oo 

P 1 = J2 np nr n ~\ 
n=2 
oo 

F' = ^2(n + 4)F n r n+3 . 

n=0 

Substituting P' and F' in the p equation one gets, 



(48) 
(49) 



Po 
6 ' 


Fi = 


_ 2pi 
15 ' 




P2 

9 ' 


F 3 - 


_ 2p 3 
21 


4P 2 pi 
105P ' 


Pi 


-H 


P2P2 


F3P1 \ 


12 ' 




18 


20 7 


2/95 


21 






27 








"( 


mpi 




2P 2 p 3 


15 


3 


7 



(50) 



This implies that po,pi,P3 can not have angular depen- 
dence whereas p n (for n > 3) have angular dependence 
provided p n -2 7^ 0. The results obtained here are consis- 
tent with the earlier work but with some changes in 
the coefficients due to 5D spacetime. Since earlier anal- 
ysis were done in QSZ spacetime, as a result, one may 
conclude that the 5D-QSZ-EGB spacetime has same local 
nakedness behavior as the QSZ spacetime. To conserve 
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space, we shall avoid replication of all other analysis be- 
ing similar to QSZ spacetime 0,0]. To further analyze 
the horizon curve, we combine Eqs. ([29)) and (|33|) giving 



t c (r)-t = 




The apparent horizon in the interior of the dust ball lies 
at R(tAH{ r ), r ) = \/F{r) — 2a. The corresponding time 

< AH( r ) is s iven b y 



t c (r) - tAn{r) 



4^2 2^5 



tan 



F -4a 



2y/2a{F-2a) 
1 



-VF -2a . (52) 



As mentioned above, at i = t c (r), we have R(t,r) = 0, 
which is the time when the matter shell r = constant 
hits the physical singularity. The singularity is at least 
locally naked if tAH > t c , and if tAH > t c , it is a black 
hole, and in case of the equality one has to compare the 
slopes. We consider the following three cases: 

a. Case I Homogeneous case, i.e., the density profile 
is homogenous and has no angular dependence, thus 



F = Fnr 4 



(53) 



Using Eq. (I5TJ1) , it can be deduced that t c < tAH (see also 
Fig.l), and hence singularity is naked. 

b. Case II Next, we assume that only F\ ^ im- 
plies that the density profile is inhomogeneous and has 
no angular dependence 



F = Fnr 4 



F ir 5 



(54) 



It is seen that t c < tAH which leads to formation of a 
naked singularity. 

c. Case III Finally, let F3 7^ which means the den- 
sity profile is inhomogeneous and has angular dependence 



F = F r 4 + F 3 r 7 



(55) 



In this case also singularity is naked because t c < tAH 
(see also Fig.l). 

Clearly, for a positive a, the central shell does not get 
trapped, and the untrapped region around the center in- 
creases with increasing a, for both homogeneous and in- 
homogeneous models. The center (r = 0) remains un- 
trapped, since for nonzero values of the Gauss-Bonnet 
parameter a > 0, above Eq. (|43l) admits no solution. In- 
terestingly, the theory demands a to be a positive num- 
ber, which forbids the apparent horizon from reaching the 
center thereby making the singularity massive and eter- 
nally visible, which is forbidden in corresponding general 
relativistic scenario. 




FIG. 1: A two-dimensional picture of the 5D-QSZ-EGB col- 
lapse showing t — t c (r) — tAn(r) versus r for the suitable 
values of parameters. The continuous curves used when den- 
sity has no angular dependence and when density has angular 
dependence dotted curves are used. 



V. EFFECT OF THE GAUSS-BONNET TERM 

It is seen here that the Gauss-Bonnet term modifies 
the time of formation of singularity, and the time lag 
between singularity formation and the apparent horizon 
formation, in contrast to the 5D dust models. Indeed, 
the time for the occurrence of the central shell focusing 
singularity for the collapse is increased in comparison to 
the 5D-QSZ case. The introduction of a Gauss-Bonnet 
term slows down the collapse process [28], which can be 
explained as follows. The contribution from the Gauss- 
Bonnet term leads, when we think in terms of a Newto- 
nian potential, to a repulsive term. The Eq. (fTo]) , for the 
marginally bound case, becomes 



R 2 



F n i? 4 



(56) 



In the a — ¥ limit we recover the 5D-QSZ solution, and 
in this limit, Eqs. (fTrJ|) and (fSrjf lead to the following 
expression for acceleration: 



Rsd- 



5D—QSZ = 



F 

R?' 



(57) 



In the 5D-QSZ-EGB spacetime the acceleration expres- 
sion (|16[) can now be expressed as 



RbD—QSZ—EGB = — 



R_ 

4a 



R_ 

4a 



l + 8a 



F 
R* 



-1/2 



(58) 



For a comparison with Newton like force equation ob- 
tained in the 5D-QSZ case $(i?) = -F/R 3 , clearly a 
positive value of a diminishes the acceleration in the col- 
lapse process. As we move out to larger shells and for 
large values of a the detailed dynamics should depend 
on higher order terms in the expansion. 

The other reason may be, for slow down of the collapse 
process in 5S-QSZ-EGB, there is relatively less mass en- 
ergy [see Eq. ([50|> ] collapsing in the 5D-QSZ-EGB space- 
time as compared to the 5D-QSZ case. This can be seen 
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from the mass function m(t,r), which is given by 

m(t,r) = R 2 (l-g ab R a R b ). 
Using Eqs. ©, 013) and $TE§ into Eq. $5§§ we get 
m(t,r) = F(r) - 2aR 4 . 



(59) 



(60) 



The quantity F(r) can be interpreted as energy due to 
the energy density e. 

Finally, for the reason mentioned above, the presence 
of the coupling constant of the Gauss-Bonnet terms a 
produces a change in the location of these horizons. Such 
a change could have a significant effect in the dynamical 
evolution of these horizons. For nonzero a the structure 
of the apparent horizon is nontrivial. In general relativ- 
ity noncentral singularity is always covered [35[ (see also 
36]). However, in the presence of the Gauss-Bonnet term 
we find that even the noncentral singularity is naked, in 
spite of being massive (F(r > 0) > 0). 



VI. DISCUSSION 

The Szekers models, depending on the sign of e, are 
subdivided in to the quasispherical, quasiplane and quasi- 
hyperbolic [l5|, [l6[ . T he g eometry of the later two is not 
really understood [HI, Il6| . On the other hand the quasi- 
spherical has been rather well investigated (Ill4l7| | , and it 
has found important applications in cosmology and grav- 
itational collapse. The QSZ metric is a dust model which 



has no Killing vector [11|, but contains LTB model as a 
spherically symmetric special case. It has been found 
that the LTB metric admits both naked singularities and 
black holes depending upon the choice of initial data. In- 
deed, both analytical [35Tj39j and numerical results [ioj] 
in dust indicate the critical behavior governing the for- 
mation of black holes or naked singularities. A similar 
situation also occurs in higher dimension the LTB models 
[3ll 434L and these results also carry over to QSZ space- 
time [3-Q. Maeda |29| and we [28| have shown that in 
spherically symmetric inhomogeneous dust collapse, the 
effect of adding a positive a does radically alter the final 
fate and leads to formation of a massive timelike singu- 
larity which is prohibited in general relativity. 

In this work, we have obtained an exact solution in 
closed form in EGB, which represents the quasispherical 
collapse of irrotational dust in 5D spacetime namely 5D- 
QSZ-EGB. The solution is nonsymmetric generalization 
of the spherically symmetric 5D-LTB-EGB solutions. It 
can be reduced to 5D-QSZ or 5D-LTB in the general rel- 
ativistic case (a — > 0). Our analysis also supports the 
earlier results [5j , deducing that under physically reason- 
able initial conditions naked singularities do develop in 
the 5D-QSZ-EGB models, which are not spherically sym- 
metric, and admit no Killing vectors. The second order 
curvature corrections changes the final fate of gravita- 
tional collapse and the nature of singularity that occurs 
in 5D general relativistic dust models. However, mild 



departure spherical collapse can not alter the standard 
picture of the structure and formation singularity, since 
the 5D-QSZ-EGB solutions discussed here are qualita- 
tively similar to the analogous spherical solutions. It is 
seen here that the Gauss-Bonnet term: (i) decelerates the 
collapse process (ii) alters the time of formation of singu- 
larities and the time lag between singularity formation, 
and (iii) modifies the apparent horizon formation and the 
location of apparent horizons. Our analysis to examine 
the nature of singularities (naked or hidden by horizon) 
is based on the comparison of tAH (time for the forma- 
tion of the apparent horizon or trapped surface) and t c 
(time for the formation of a central singularity). In QSZ 
the singularity can be directionally naked Q. Hence, it 
would be necessary to investigate in further details the 
final fate of the inhomogeneous dust collapse in the EGB 
theory in order to bring out explicitly the difference in 
global nakedness when we depart from spherically sym- 
metric [42l | . 

The conjecture that such a singularity from a regular 
initial surface must always be hidden behind an event 
horizon, called cosmic censorship conjecture (CCC) was 
proposed by Penrose [4l[. The CCC forbids the exis- 
tence of naked singularities. Despite almost 30 years of 
effort, we are far from a general proof of CCC (for re- 
cent reviews and references, see [431] ). But, significant 
progress has been made in trying to find counter exam- 
ples to CCC. Our analysis, as in the 5D-LTB-EGB [28|. 
shows that there exists a regular initial data which leads 
to a naked singularity and hence in our nonspherical case 
also the CCC is violated. The usefulness of these mod- 
els is that they do offer an opportunity to explore the 
properties of singular spacetime. The investigation of a 
mild departure from standard spherical symmetry mod- 
els may be valuable in attempts to put CCC in concrete 
mathematical form. Finally, studying such models which 
lacks symmetry is important, so that one can check which 
properties of gravitational collapse are preserved. This 
may also helps us to bring out some universal features in 
the theory of gravitational collapse (42| ■ We have shown 
here that there exist regular initial data which leads to 
a naked singularity violating CCC. However, this may 
not be a serious threat to CCC because of the following 
two reasons viz. (i) The matter considered here is dust 
which is only an effective, macroscopic approximation to 
a fundamental description of matter 44] . (ii) In general 
relativity the energy-momentum tensor given by Eq. (|8[) 
satisfies the weak energy condition. However, this may 
not be true in EGB because the Gauss-Bonnet term itself 
violates the energy condition. 
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